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I. INTRODUCTION 

This report summarizes the results of a set of observations of the 
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atomic oxygen 01 15,867 K 6300 8 V- D thermospheric emission to determine 
the temperature of the emitting species using a Fabry-Perot interferometer 
at the University of Alaska Geophysical Institute (64.86° latitude, 
-147.85° longitude) during March 1984. Spectral profiles obtained from the 
interferometer are used to determine the Doppler temperature by means of 
the technique reported by Hays and Roble [1] and Roble [2]. 

This general concept of using Doppler widths of airglow lines to find 
temperatures has been used by many investigators such as Wark [3], Nllson 
and Shepherd [4], Turgeon et al. 15], Zwick and Shepherd 16], Hernandez 
[7], Hernandez and Roble [8-11], and Smith et al. [12]. Analytical 
descriptions of Fabry-Perot spectrometers have been presented by Born and 
Wolfe [13] and Hernandez [14]. A review of temperatures and winds measured 
by Fabry-Perot spectrometry was done by Hernandez [7]. 

II. INSTRUMENTATION 

The Fabry-Perot interferometer was located on the top floor of the 
Geophysical Institute, Fairbanks, Alaska, allowing observations through a 
variable geometry periscope system and a Plexiglass dome. The interfero- 
meter is described by Roble [2] and Sivjee et al. [15] and is shown in 
Figures 1 and 3. Spectral profiles were obtained by pressure scanning of 
the etalon in steps. A stepper motor controlled the pistor (shown in 
Figure 2), and position data were provided by an A/D encoder. During a 
pressure scan, the refractive index in the etalon changed. This caused the 
light in the center of the interferometer pattern to change intensity, 
creating an intensity /pressure fringe. This light from the interferometer 
was reflected through a 1/8-inch aperture, then filtered, and finally 
detected by a photomultiplier tube. Simultaneously, a photometer recorded 
the Intensity of the 6300 R lines in the same direction as the 
interferometer observation. Also, the etalon temperature was monitored. 
Pressure scanning, mirror positioning, and integration times were 
controlled by computer. Times, pressures, observation directions, 




temperatures, and counts were automatically stored in the format shown in 
Appendix F. 

An He-Ne (6328 X) frequency-stabilized Tropel laser was used for 
calibration. To prevent photomultiplier tube saturation, the beam was 
incident on a ground glass sphere above the interferometer. The glowing 
sphere simulated the sky, so pressure scanning produced fringes without 
changing the experimental setup from the sky observation configuration. 
If, during calibration, the fringes were found to be asymmetric, the 
aperture or, in extreme cases, the etalon plates were adjusted. Neither of 
these adjustments were required between the observations included in this 
report, although such adjustments were required between earlier trial runs. 


III. OBSERVATIONS 

Observations were made at an inclination of 30° from the horizontal in 
the north, south, east, and west magnetic directions. Observations were 
also made at local zenith. All observations were made between sunset and 
sunrise from March 6 to March 16, 1984. At the beginning of this period 
the moon was only 9% full; but by March 16, it was 97% full. 

Most nights were clear, and slight to moderate auroral activity was 
observed in the zenith and to the north. On clear nights, diffuse aurorae 
were usually seen to the north. Data were also taken on cloudy and foggy 
nights, even though the observed light was scattered from all sky sources 
prior to detection. If there is mass average velocity with the emitting 
atoms, the temperature measurement is still correct. Winds could cause an 
error if they were unusually strong under these cloudy conditions. 

In addition, the raw data were monitored during all runs. For 
example, an X-Y plotter was used to record the intensity of the 
interferometer image as a function of pressure to monitor the fringes and 
to provide a quick look at the raw data. After each night of observations, 
the records stored in the computer were transferred to disks. The 
parameters that were monitored and the formatting of this data are 
discussed in the following section. 
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IV. DATA REDUCTION 


The computer program for data analysis included: reading of data from 

disk, temperature and pressure compensations, normalization, Fourier 
transform smoothing, deconvolution, and least squares fitting to 
theoretical profiles. A listing of this program is included in Appendix D 
along with a sample output. 


A. Data Format, Reading Data 

Each pressure scan of 64 steps included one free spectral range or 
fringe, which covered only 34 steps. The other 30 steps corresponded to 
other fringes and were, therefore, ignored. These fringes occurred as 
variations in photomultiplier counts as functions of pressure. Data at 
each step included etalon pressure and temperature, photometer reading, and 
photomultiplier reading. Times were recorded at the beginning of each 
pressure scan. All these data were stored in files such as the one shown 
in Appendix F. The first step in the data analysis program was to read the 
data corresponding to one fringe. 

B. Data Modification 

First, the reading of the etalon pressure was adjusted to compensate 
fcr variations of etalon temperature. Then the photomultiplier counts 
(i.e., intensity) were adjusted to account for temporal intensity 
variations of the sky by dividing by the photometer reading, which was 
recorded in kilorayleighs. Finally, the background signal was subtracted. 

C. Fourier Cosine Transform Smoothing 

The Fabry-Perot interferometer produces a series of concentric rings. 
The center of this image is monitored by a photomultiplier tube and as the 
etalon pressure is changed, the measured intensity changes. A series of 34 
different pressure readings were taken in each scan. The resulting 
observed profile is a convolution of the actual sky profile and the 
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instrument function, with noise superimposed. To extract this noise a 
Fourier technique was used. 

A Fourier cosine transform is fit to the data from one fringe. When 
the first six coefficients are used to reconstruct the curve, the result is 
a smoothed version of the data as shown in the equations below. If not 
enough coefficients are taken, the data are not accurately represented. On 
the other hand, if too many coefficients are taken, the resulting curve 
fits statistical noise as well as the source variation. 

For a theoretical emission profile which is perfectly smooth, the 
coefficients decrease as the wavenumber (i.e., m in equation (1)) 
increases, as shown by Hays and Roble [1]. For our data, the coefficients 
decrease until m is 6, whereupon it begins to increase, as shown in Figure 
4. Therefore, to represent the actual data with a curve that is as accurate 
as possible with minimum noise, the optimum number of coefficients is 6. 
Each of the three curves shown in Figure 4 are from different pressure 
scans, all taken successively; therefore, each curve corresponds to a 
different fringe. 

The smoothed profiles of the observed data are given by: 
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Y = E Y cos (mx) , 
m=l m 


( 1 ) 


where the Y 1 s are the Fourier coefficients from fits to the observed data, 
m 

These coefficients are found using: 


Y = 1/ttv/y 2 + Y T , 
m cm sm 

where 

- l] £ [2 sin (c2 P ) (sin c2 ^ ) 2 
cm i=l m i 2. 

+ cos (c2 P^) sin (c2 DPS)] 
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where 
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C = counts at ith step 
= pressure at ith step 


2tt x m 
SRS 


SRS = free spectral range, pressure units 
DPS = P t - P t-1 ; 

and similarly, 


Y 

sm 


34 

\h 


C 

m 


12 cos (c2 P t ) (sin c2 ^) 2 


- sin (c2 P^) sin (c2 DPS)] . 

D. Deconvolution 

The smoothed profile is assumed to be a convolution of the theoretical 
Doppler profile of the sky emission (D) and the instrument profile (L) 
(obtained by observing through the instrument a diffuse illumination 
provided by an He-Ne laser); i.e.. 


Y = D * L . 
m m m 

To deconvolve the laser profile from the observed data, the 
coefficients from the observed data were divided by the laser coefficients 
using the convolution theorem (e.g., see Reference 16); i.e., 


m 


Y /L , 
tn m 
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where Che laser coefficients are found by running a slightly modified 
version of the program (In Appendix D) independently for the laser fringe 
coefficients only. 


E. Least Squares Fit 

After fitting the data to a smoothed profile using the Fourier 
coefficients and deconvolving the instrument function, the profile is then 
compared to a series of theoretical Gaussian emission profiles. These 
theoretical profiles are discussed in Appendix A. They have a known 
Doppler width for any given temperature. These theoretical profiles are 
convolved with a Gaussian instrument function obtained by calibrating with 
an He-Ne laser (as in equation (2)). The theoretical curve that gives the 
least square error with the smoothed data was used to obtain the 
temperature. 

To ensure that the deconvolved profiles are Gaussian, they are plotted 
with theoretical Gaussians as shown in Figure 5. Non-Gaussian results were 
. ~’t included in temperature plots. Also, the half widths of all 
deconvolved profiles have been compared to those in Table A-l, Appendix A 
to ensure that no major errors occurred in the deconvolving calculation. 


V. RESULTS AND DISCUSSION 


After converting the fringes to temperatures, the temperatures were 
plotted as a function of time for each night of observation, as shown in 
Figures 6a-g. Each direction of observation: zenith, north, south, east, 

and west (all magnetic) were included. The times given are in Universal 
Time (UT) , where local time was 9 hours behind UT. 

It Is apparent from the results in Figures 6a-g that the data points 
are scattered, but are typically within a band that is ±100 K about the 
mean. The data for March 7 (Figure 6a) are highly scattered, but on March 
10, the temperature is seen to be about 800±100 K. The next three nights 
give temperatures of 750, 750, and 800 K. All the above data were taken on 
clear nights. March 14 was cloudy and the apparent temperature was found 
to be 750 K. March 15 was foggy and the apparent temperature was found to 
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be 700 K. On March 16, only four zenith points were taken and three of 
these were at about 680 K. 

Contributions to the spread in the data may be from both the 
instrument and actual variations in 01 emission intensity. It is difficult 
to determine the fractional contributions of these components, although it 
has been considered by Hays and Roble [1]. 

A. MS1S 


Also plotted in Figures 6a-g are the temperature profiles derived from 
the MSTS (Mass Spectrometer Incoherent Scatter) model of the neutral 
atmosphere [17-18]. Parameters utilized in this model include local time, 
altitude, location, 3-month average Fjq ^ flux, dally Fjq ^ flux, and daily 
Ap value (see Appendix B) . An example of the MSTS output is given in 
Appendix C. 

The purpose of running MS1S was to deduce the predicted model 
thermospheric temperature for the geomagnetic conditions at the observing 
location. Comparisons of the MSIS temperatures with our results show that 
ours are in a reasonable range. 


TABLE 1. PLOT SUMMARY 


Date (UT) 

MSIS Height 
(km) 

Temp. 

(K) 

Sky 

March 7 

high dispersed 

clear 

Maich 10 

180 

800±100 

clear 

March 1 1 

170 

750+100 

clear 

March 12 

170 

750+100 

clear 

March 13 

170 

800+100 

clear 

March 14 

170 

750+100 

cloud 

March 15 

150 

700+100 

fog 

March 16 

minimal data 

680 

clear 


B. Shadow Height 


In order to determine whether the observed volumes were sunlit, the 
Earth's "shadow height" was calculated. Heights for zenith, north, south, 
east, and west (magnetic) observations were c&lculated. It was determined 
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that the observed volumes In some of the observations (taken shortly after 
sundown and shortly before sunrise) were sunlit. Plotted in Figures 6a-g 
are heights of the Earth’s shadow. This shadow height is the distance from 
the point where the line of observation intersects the surface between 
sunlight and shadow to the ground directly below. (See Appendix E for a 
detailed explanation of how shadow height is determined.) 

Each curve corresponds to a given direction of observation (e.g., 
east). Where the curve intersects an MSIS altitude curve, the 
corresponding time indicates when the shadow height was at the altitude of 
the MSIS curve. For example, in Figure 6a, the vertical curvr shown 
corresponds to observation toward the east. This curve intersects the MSIS 
170 km curve at 1400 UT, indicating the shadow height at that time. As 
shown, the shadow height passed through lower altitudes at later times. 
The overall curve indicates sunrise, because observation points to the left 
of the curve are in darkness, while observation points to the right, i.e., 
later in time, are sunlit. In the following figures, different directions 
of observation are included. N represents north, while W and E represent 
west and east, respectively, Also shown in the figures are the curves for 
sunset. Of all directions observed, points to the west were the last to 
remain sunlit on any given night, and points to the east were the first to 
become sunlit. 

For sunrise observations, the lowest shadow heights were to the east 
It can be seen that on March 11 and 13, observations to the east, an hour 
after the emission region passed into sunlight, indicate pre-dawn 
temperature enhancement. After sunset, the lowest shadow heights were to 
the west, but sunlit dusk observations were made only on March 14, which 
was a cloudy night. 

Shadow height calculations were based on Chamberlain’s work Tl9l. 
Figure 7 shows the resulting shadow heights as a function of time for all 
directions observed. Programs In Appendix E were used to obtain these 
results. 


C. Volume Emission Rate Profiles Previously Measured 

Altitude profiles of the 6300 X volume emission rate have been 
reported by various investigators using data from the Atmosphere Explorer 
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program. Selected results are noted here for comparison with the results 
presented in the previous section. For example, Abreau et al. [20] found 
altitude profiles of the volume emission rate using a photometer onboard 
the AE-E satellite. From November 1980 to February 1981, averaged data 
gave a peak rate ax. altitudes of 250-280 km. 

Hays et al. [21] reported altitude profiles of the 6300 A volume 
emission rate for various solar zenith angles during evening twilight. 
Peak emissions were in the 220-250 km altitude region. Their theoretical 
models using O 2 photodissociation and phot' lectron impact indicated that 
below 200 km, the dominant mechanism of excitation was photodissociation. 

Mid-latitude orbit data for summer and winter were examined by Torr et 
al. [22]. Theoretical and measured peak emission altitudes were 
approximately 200 km in winter and 200 km in summer. 

Hernandez and Roble [10] obtained temperatures from the 6300 X 
emission as a function of time. Their temperatures were found to be 
greater than MS1S model predictions for an altitude of 250 km during summer 
and equinox nights. Their data for March 1976 (the same month as our 
data) were scattered to about the same degree as our data, and their data 
were centered on 800 K in agreement with our results. Further results were 
later reported by Hernandez and Roble [11]. They found that the 6300 %. 
line peak emission rate occurred below the F 2 peak at 250 km, at about 100 
R, decreasing throughout the night. 

Recently, Sipler et al. [23] also measured the neutral F-region 
temperature using a Fabry-Perot Interferometer. During geomagnetically 
quiet nigh' : from 1975 to 1979, equinox solar minimum average temperatures 
were found to be about 750 K. 

Emission rate profiles for 6300 dayglow were reported by Killeen et 

—3 —1 

al. [24], They reported that a typical profile had a peak of 180 cm s 
at about 210 km altitude. 

The above results indicate a peak volume emission rate at altitudes of 
200-300 km. As discussed earlier, these results are in rough agreement 
with the present observations, given the spread in our temperatures and the 
combined uncertainties of MSIS temperatures, as well as our lack of 
statistical data. Assuming the glow intensity is originating from 200 km, 
our data suggest the MSIS model is predicting slightly higher temperatures 
than observed. 
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VI. FUTURE WORK 


Many improvements could be made to the apparatus. With CCD arrays now 
readily available, it is possible to image the Fabry-Perot fringes. 
Imaging of fringes has been done by Sivjee et al. (15] and Rees and 
Greenaway [25]. The latter investigators developed a Doppler imaging 
system (DIS) which used a Fabry-Perot interferometer and a 120° field of 
view all-sky camera. The multiring image contains spectral and spatial 
information. By imaging the fringe from the interferometer (of this 
investigation), all tne information could be digitized and analyzed by 
computer. By using the entire image of the fringe rather than just the 
central fringe, a much stronger signal is obtained. This system would 
provide a convenient means of monitoring winds using the Doppler shift. 
Yet another advantage of this system would be that fringes would be made 
much more quickly than with pressure scanning. Further, integration times 
could be easily varied to an optimum duration. This system could be 
further enhanced by using a pyramid with mirror surfaces rather than a 
scanning mirror. This would allow monitoring of all four directions 
simultaneously. With both of the above modifications, all moving pares 
would be eliminated, greatly simplifying the entire system and reducing its 
size. 

An extension of the pyramid system would be to use an all-sky lens so 
that all directions could be monitored simultaneously. This lens could be 
coupled to an interference filter, an image intensifier, and finally to a 
CCD array via a fiber optic plug. In this arrangement, the apparatus could 
be considerably reduced in size compared with the current system. However, 
as size decreases, integration time increases. An image intensifier could 
be used to reduce integration time. 

Cooling of the array would be imperative to reduce noise accumulation. 
A Peltier electric cooler, for example, would be convenient for this task. 


VII. CONCLUSION 

The temperatures of the 6300 line on the clear nights of March 10, 11, 
12, and 13 were 800, 750, 750, and 800 K, respectively, with a spread of 
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±100 K. Despite the large spread in the data, the temperatures and the 
spread in the temperatures were consistent from day to day, and consistent 
with previous observations. 
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Figure 1. Fabry-Perot Interferometer control schematic 









Figure 3. Fabry-Perot Interferometer controller and data acquisition and display system schematic. 









Figure 5, Comparison of theoretical Gaussian profiles to smoothed data 
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Figure 6e. March 13 
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Figure 7. 


Shadow height, as a function of Local Time in different look directions 
for days encompassing the period of the observations. 
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APPENDIX A. THEORETICAL GAUSSIAN PROFILES 


A Doppler shift of wavelength results from the motion #f a radiating 
particle towards or away from an observer. In gas or plasmas, the random 
thermal motions of all particles lead to a Maxwellian velocity 
distribution. This results in a Gaussian distribution in observed 
frequency due to the Doppler shifting. Hence, the spectral profile has a 
Doppler broadened component that is a function of temperature. From Griem 
[26], the Doppler shape is given by 


I ( A A ) = 


m 2 h 

( 

c 

exp < 

2TrktA^ 

o 

( 


-M 


2kyA 


t(aa) 


2 


* 


(A-l) 


o _7 g 

where A » 6300 A = 6.3 x 10 m, c = 2.99792 x 10 m/s, M = mass of 
o — 26 

emitting species (0) (= 2.6776 x 10 kg) , k = Boltzman's constant (= 1.38 

oo 

x 10 J/K), and T = temperature in K. 


Normalized intensity as a function of T and AX is found from: 

,26 (AX ) 3 


Normalized Intensity = exp 


-2.19682568 x 10 


The above equation is verified as shown in Appendix G, by comparing it with 
Wark's [3] version of this equation. 


theoretical Gaussian Profiles 


To convert the form of this equation to one in terms of etalon 
pressure instead of wavelength, X, the conversion 


A> 


1.98 x 10 -11 mFSR 

— — — : — — x A pressure units 

74.09 pressure units FSR v 


(A-2) 


is substituted into equation (A-l) to obtain: 
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Normalized Intensity ® exp 


-15.689416 


( A pressure units) z 
T 


where FSR is Free Spectral Range and T is temperature in K. 


The factor (1.98 x 10 ^ m FSR) is a fixed instrument parameter. The 

FSR observed with the calibration was 75.53 pressure units, and this value 

is used throughout the program and in the calculations except in the above 

step. Here, a slightly smaller FSR was used to compensate for the lact 
o 

that the 6328 A laser was used t^r calibration, while the observations wsre 
o 

of 6300 A. Since the FSR is proportional to the square of the incident 

wavelength, the 6328 He-Ne laser which gave a FSk rf 75.53 would, at 6300 
o 

A, have given an FSR of: 

7 5.53 (||g)* - 74.86 

So the FSR of 75.53 was 75.53 - 74.86 = 0.67 pressure units too high. To 
compensate for this, the FSR used in the above step (equation (A-2)) is 
approximately 75 pressure units. 

To input A pressure units in radians x 100 (where 2ir rad of 
trans; >rmed data corresponds to 1 FSR of data), the conversion for each set 
of data 

75.53 pr e ssure units _ 12.02097 , _ 

2tt rad x 100 * 100 


is substituted into equation (A-2) to obtain: 


Normalized Intensity - exp 


0.2267 * ad)2 | > 


where the numerical factor 0.2267 is labeled "CO'* in the program. The 
reason for scaling the pressure units by the factor 100 is to make a unit 
change in input (equation (A-3)) small enough to allow many (628 'v 100 x 2 
x if ) steps per fringe. This facilitated programming of the graphics. 

To check the results of the program, the HWHMs of the laser and 
theoretical Gaussians can be easily convolved. The laser width of 3.5 
pressure units was obtained from a pressure and temperature compensated 
laser calibration scan. To find the HWHM of the theoretical Gaussians 
start with equation (A-2). 
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Normalized Intensity = exp |-15.689 k re, - gf j - 

Let normalized intensity » 0.5 (i.e., HWHM) , then solve for pressure units 
to obtain: 


HWHM in pressure units 


T In 0.5 
-15.689 


Convolving this with the laser calibration scan width provides a quick look 
determination of the temperature from the observed width. Examples are 
shown in Table A-l . 


TABLE A-l. QUICK-LOOK DETERMINATION OF TEMPERATURE FROM OBSERVED WIDTH 


Temperature 

(K) 

Theoretical Gaussian HWHM 
(pressure units=psi) 

Convolution of Theoretical 
Gaussian and Laser/llWHM 2 + 3.5 Z 
(pressure units**psi) 

1000 

6.646 

7.51 

1100 

6.97 

7.80 

1200 

7.28 

8.078 

1300 

7.57 

8.3477 

1400 

7.86 

8.608 

1500 

8.14 

8.8611 

1600 

8.407 

9.107 

1700 

8.666 

9.346 

1800 

8.917 

9.5799 

1900 

9.162 

9.8077 

2000 

9.400 

10.030 
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APPENDIX B. SOLAR GEOPHYSICAL DATA 
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APPENDIX C. MSIS EXAMPLE OUTPUT 


MARCH 7 UT 
84067 

INPUT UNIVERSAL TIME IN SEC 
14440 (0400 UT) 

INPUT MIN, MAX ALTITUDES AND ALT. STEP IN KM 
130,290 20 

INPUT GEODETIC LAT., EAST LONG. IN DEG 

64.9.212.2 

INPUT LOCAL APPARENT SOLAR TIME IN HRS 

18.2 

INPUT 3-MO. AVE. OF F.g ? FLUX 

112 

INPUT DAILY F. n - FLUX FROM PREVIOUS DAY 

109.5 

INPUT DAILY AP VALUE 

26 

DATE UT(SEC) LATITUDE LONGITUDE LOCAL TIME <F _ > F in7 AP 
84067 14440 64.9 212.2 18.2 112.0 109.5 26 


ALT [HE] [0] [N2] [02] [AR] [H] T 


130. 1 . 7502E+07 3.2964E+10 1.1277E+11 2.0486E+10 4.8818E+08 1.4128E+06 524. 
150. 1.2866E+07 1.3711E+10 3.0166E+10 4.7443E.09 8.6045E+07 5.0339E+05 692. 
170. 1.0393E+07 7.1768E+09 1.1209E+10 1.5684E+09 2.2538E+07 2.5966E+05 806. 
190. 8.8380E+06 4.2385E+09 4.8880E+09 6.1621E+08 7.1885E+06 1.7916E+05 883. 
210. 7.7383E+06 2.6843E+09 2.3331E+09 2.6699E+08 2.5662E+06 1.4700E+05 936. 
230. 6.8957E+06 1.7755E+09 1.1773E+09 1.2287E+08 9.8243E+05 1 3188E+05 972. 
250. 6. 2138E+06 1.2081E+09 6.1598E+08 5.8823E+07 3.9375E+05 1.2358E+05 996. 
270. 5.6406E+06 8.3788E+08 3.3030E+08 2.8927E+07 1.6285E+05 1.1829E+05 1013. 
290. 5. 1460E+06 5.8885E+08 1.8022E+08 1.4498E+07 6.8877E+04 1.1441E+05 1025. 
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APPENDIX D 


PROGRAM AND EXAMPLE OUTPUT 


•TYPE G.G 

C PROG TO GENERATE THEORETICAL GAUSSIAN PROFILES 
C AND PLOT OBSERVED PROFILES AND FIND TEMPERATURE 
C Following is the starting point ( sp ) value used to 

C establish the place in the data to start reading* 

SP-1668. 

C Following is a table that appears in the report? and is not essential 

C to this program* It is included here for reference only* 

C GENERATE A TABLE TO FIND APPROX TEMP FROM HUHM OF RAW DATA 

CC TE*700. 

CC WRITE( 5f # y ******* 

C X IS TEMP? XX IS HUHM PSI? XXX IS HWHM CONVOLUTED WITH 3*5 PSI LASER 
CC URITE<5?*>'TEMP? HWHM THEORETICAL HWHM CONVOLUTED ' 

C3 X*( ( TE*< ALOG< *5 ) > >/-15*689 >#**5 

C XX=X*3l5*/27*67 

C XXX=< ( X*X >+< 3*5*3 ♦ 5 > >*#*5 

C WRITE< 5?*)TE?X?XX?XXX 

C TE=TEf 100* 

C IF (TE *LE* 1700) GO TO 3 

C ******* ***************************************************************** 

C The program starts here* YSM? YSSM? and YSCM are the founer cofficients 

DIMENSION YSM( 15)?YSCM< 15)?YSSM( 15) 

REAL F< 10 ) 

C The following step is reouired for graphing only and is not essential 

CALL IN1TT( 30 ) 

C RETREIVE OBSERVED DATA 

C RETREIVE OBSERVED TEMP(TL)? PRESSURE< PL )? COUNTS(CL)? 

C KILORAYLEIGHS< KL )? NP IS NEW PRESSURE? NC IS NEW * OF COUNTS 

REAL TL( 34 ) 

REAL PL( 34 ) 

REAL CL( 34 ) 

REAL KL( 34 ) 

REAL NP< 34 ) 

REAL NC( 34 ) 

0PEN< UNIT*3?NAME*' FW2JFP16.DAT' ? READONLY ?TYPE=' OLD' ) 

c read: temp pres counts krays 

C READ BLANKS UNTIL START OF DESIRED DATA********************* 

ZW«1 

URITE< 5?*)'ZU ZS' 

187 READ( 3?*?ERR-187 ) 

C WRITE( 5?*)ZW?ZS 

ZW»ZU+1 

C MAX ZW VALUE (SP) BELOW DETERMINES STARTING PT TO READ DATA 
C For convenience? the following step was moved to near begining of prog* 

C SP=1185. 

IF (ZW *LE* SP) GO TO 187 

ZY=1 

ZY=ZY+1 

C C ************************************************************* 
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C BlanKs have now been read to begining of data* Start reading data now* 
ZZ-1 

189 READC 3*150 >TLC ZZ >* PLC ZZ >* CLC ZZ )* KL< ZZ ) 

C PRESSURE AND INTENSITY COMPENSATION **♦♦*#♦*♦ 

NP< ZZ )»PL( ZZ )*< 298. /< 273.4-TLC ZZ ) ) ) 

NC( ZZ >-CL( ZZ >/KL< ZZ ) 

C **********♦*♦#**♦**♦♦♦♦♦♦♦♦♦**♦♦*♦♦♦♦*♦♦♦♦♦♦ 

ZZ*ZZ+1 

C The following step tells computer to read 34 lines of data (the fringe). 
IF <ZZ *LE* 34) GO TO 189 

150 FORMAT < 9XpF6*3*F7*2#13X»F6*0f4XfF5*3 ) 

C CALCULATE BACKGROUND INTENSITY ZB-AVERAGE OF 1ST AND LAST PTS 
C OF NC AFTER COMPENSATIONS* BEFORE SUBTRACTION OF BG AND NORMALIZATION 

ZB*< NC( 1 HNCC 34 ) )/2 
WRITEC 5**)' ttttttttttt' 

WRITE( 5*# )' BACKGROUND INTENSITY IS' 

WRITE( 5p# )ZB 

Following step not required* so ignore* 

WRITEC 5*# >' NP( 34 )-NP< 1 )»' 

Following is specif icstion of free spectral range in pressure units 

SRS*75.S3 

WRITEC 5r* )SRS 

NORMALIZE NUMBER OF COUNTS AND SUBTRACT BACKGROUND 

NORMALIZE BY DIVIDING BY 11TH ELEMENT OF NC ( 18TH FOR 34ELEMENTS ASBELOW ) 
WRITEC 5*# ) ' NC BEFORE NORMALIZING AND AFTER SUBTRACTION OF BG' 

ZA*1 

149 WRITEC 5**> NCC ZA )-ZB 
ZA*ZA+1 

IF ( ZA * LE * 34 ) GO TO 149 

NORMALIZE TO NN»NC< 11 )-ZB®MlDDLE NO. OF COUNTS - Z BACKGROUND 
WX*NC( 18 )-ZB 
WRITEC 5**)'WX IS' 

WRITEC 5*# >WX 

The next 4 steps are used as checKs only and may be ignored* 

NN-NCC 18 )-ZB 
NN°WX 

CC WRITEC 5*#)' WX IS' 

CC WRITEC 5*#)WX 

C The next 4 lines are the normalization described 15 lines above 
ZU»1 

152 NCC ZW >«( NCC ZW >-ZB )/WX 
ZW=ZW+1 

IF (ZW *LE* 34- ) GO TO 152 

WRITEC 5*# )' PRESSURE* COUNTS* KILORAYLEIGHS* IN RAW DATA FORM' 

WRITEC 5*# )' LAST 2 COLUMNS ARE NEW PRESSURE AND NEW COUNTS' 

WRITEC 5 **>' PL CL KL NP NC' 

ZA*1 

151 WRITEC 5**>PL( ZA),CLC ZA )*KL( ZA )*NP< ZA)*NCC ZA ) 

7 As 7A11 

IF ( ZA .LE. 34 ) GO TO 151 

URITE(5.#)' THIS TEXT CAUSES LAST LINE TO PRINT' 

C 100++++++++-H-+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


(JUU O O CJ 


k < s 






C Newt two lines are free spectral ranse in pressure units 
URITE(5.*>'SRS' 

WRITE< 5»S )3R8 

Newt line applies on la to sraphina on terninal 
CALL FINITTt 0.760 ) 

NOW FIND COEFFS OF OBSERVED DATA ttttttttttf ttttttttTtf Tttttttttt 
PI-3.1415926535 
URITE<5»*>' ' 

The followins is used if the Delta Pressure Step is desired output 
URITE<5»«>' DPS NP< J )' 

MAX M IS NO. OF FOURIER COEFFS 
DO 106 M-1,10 
C FSR IN PRESSURE UNITS 
C SRS-75.53 
SUNC-0. 

SUNS-0. 

C2-6.28318IKM/SRS 

DO 107 J-l »34 

IF (J .NE. 1) 60 TO 104 

C Followins line calculates the difference in pressure between 1st 2 steps 
DPS-NP< J+l >-NP< J ) 

C DPS»< SRS-( NP( 20 >-NP< 1 ) ) )/2 
60 TO 105 

104 CONTINUE 

C Followins line calculates the difference in pressure between steps 
DPS«NP( J >-NP< J-l ) 

C IF (J .NE. 20) CO TO 105 
C DPS-( SRS-< NP( 20 >-NP< 1 ) ) )/2 

105 CONTINUE 

IF ( M .NE. 1 > 60 TO 116 

C The *ollowins step is used if Delta Pressure Step is desired output. 

C URITE( 5.4 )DPS»NP( J ) 

116 CONTINUE 

C The newt two steps calculate the Fourier cos coeffs. and sin coeffs 
SUNC»SUNC+< NC( J )/M >*< 2.*SIN< C2*NP< J ) )*( < SIN< C24DPS/2. ) )**2 
* HCOSt C2#NP< J > >*SIN< C24DPS > ) 

SUNS-SUNSH -NC( J >/M >*< 2.*C0S< C2*NP( J ) >*< ( SIN< C24DPS/2. > >**2 
* )-SIN< C2#NP< J > >*SIN< C24DPS ) > 

107 CONTINUE 

C The newt line calulates the Fourier coeffs of the frinae 
YSM< M >»SQRT( SUNS*SUNS+SUNC*SUNC )/PI 
YSCM< M >»SUNC/PI 
YSSM< M J-SUNS/PI 

106 CONTINUE 

C Followins are the Fourier coeffs 
UIRITE( 5.4 )' YSM' 

URITE<5.4)YSM< 1 ) 

WRITE(5»4)YSM<2> 

WRITE( 5.4>YSM< 3) 

WRITE( 5.4>YSM< 4) 

URITEC5»4>YSM<5> 

URITE(5»4)YSM<6> 


34 





c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

900 


C910 

c 

c 


c 


c 

c 

c 

c 

C920 


2 ? 

*>00# 


p An 


URITE< 5t*)YSM<7> 

URITE< 5»* )YSM< 8) 

URITEC 5f#)YSM< 9) 

WRITE< 5r*)YSM< 10) 

WRITE< 5*# )YSM< 11 ) 

URITE< 5» * )YSM( 12 ) 

URITE< 5t*)YSM< 13) 

WRITEC 5,* >YSM< 14) 

WRITE< 5*# )YSM< 15) 

BELOU IRRELEVANT EXCEPT FOR LASER COEFFS**#**#»**###**#**#*##**#* 
UJRITE< 5 *# V NOW PLOT OBSERVED DATA USING FOURIER CQEFFS' 

CALL INITT( 30 ) 

CALL DWINDO< 0 . » 629 ♦ r-l*>l. ) 

CALL TUINDO< 50*900*50*800) 

CALL MOVEA( 0 • *0* ) 


IS 

Quality 


X-.01 

FIND FIRST VALUE FW AND CENTRAL CU TO N PLOT ONLY 
FU*YSM< 1 )*COS< .01/100 >4YSW 2 >*COSI 2*. 01/100 ) 

#+YSM< 3 )#COS< 3*. 01/100 )+YSM< 4 >*COS< 4*. 01/100 ) 

#+YSM< 3 >*COS< 5*. 01/100 )+YSM< 6 >*COS< 6*. 01/100 ) 

#+YSM< 7 )#COS< 7*. 01/100 >4YSM< 8 )*COS< 8*. 01/100) 

CU»YSM< 1 )*COS< PI >+YSM< 2 )«COS< 2*PI >+YSM< 3 )«COS( 3#PI ) 

#+YSM< 4 >*COS< 4*PI >+YSH< 5 >*COS< 3*PI )+YSM( 6 >*COS< 6#PI ) 

#+YSW 7 >*COS< 7*PI )4YSH( 8 >*COS< 8*PI >+YSM< 9 )*COS< 9*PI ) 

Y*< ( YSM( 1 >*COS( X/100 )+YSM( 2 >*COS< 24X/100 )4YSM( 3 >*COS< 36X/100 ) 

* + YSM< 4 )*COS( 4AX/100 >+YSM< 5 >*COS< 5*X/100 )+YSM( & )*COS( 6*X/100 ) 
#+YSM< 7 >*COS< 7*X/100 )+YSM( 8 )*COS( 8*X/100 )+YSH( 9 )#COS< 9*X/100 ) 
#+YSM< 10 >*COS< 10KX/100 >+YSM< 11 )#COS( 11#X/100 ) 

#+YSM< 12 )*COS< 12*X/100 >+YSM< 13 >*COS< 13*X/100 ) 

# + YSM< 14 >*COS< 14*X/100 >+YSM< 15 >*COS< 15*X/100 ) ) ) 


CALL HRAWA(XfY) 


X=X4-1 

IF (X .LE. 629) CO TO 900 

NOW PLOT LASER PROFILE WITH FOURIER COEFFS 
X=.01 

REAL L( 10 ) 

LASER COEFFS FOR FSR OF 75.7332 PSI 

L< 1 )= . 2254809 

L< 2 >-.1613598 

L( 3 >*.1071364 

L< 4 >=.094267033 

L( 5 >=.074299119 

L( 6 >=.055807292 

L< 7 >=.029683163 

L( 8 )= . 02P19389 

L< 9 >*.018680153 

L< 10 >*.018559434 

200F++++++++++++++++++++4«-++++:-+++++4+++++++++++++++4+++++++4++4+4+++++ 
LASER FIRST VALUE ( LF ). LASER CENTRAL VALUE < LC > TO N PLOT ONLY 
FF=. 7 633933 
CC=-0. 1496794 

Y=( < L< 1 )#COS< X/100 >+L< 2 )*COS< 2#X/100 >+L( 3 >#COS< 3*X/100 ) 
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920 Y-LC 1 >*C0SC X/100 )+L< 2 )#C0S< 2*X/100 )4LC 3 >*CQS< 3*X/100 ) 

#+L< 4 )*C0SC 4*X/100 HLC 5 )*C0SC 5*X/100 ) 

*4L< 6 )*C0SC 6*X/100 )*LC 7 )*C0SC 7#X/100 ) 

*4LC 8 >*C0S< 8*X/100 HLC 9 )*C0SC 9#X/100 ) 

#4*L< 10 )*CQS< 1C X/100) 

CCC CALL PRAWAC Xf Y ) 

X-X+l 

CCC IF <X . LE. 629) GO TO 920 

C *#**#**##•**##*#♦*#######*«#**##*##«##«*##♦####***#*##*♦*•♦*#•*# 

C 10 Fourier coeffe were found but anu number up to 10 mau be used 

C Now select the number of Fourier coeffe to be used: 

C YSK, 3)«0. 

C YS.i( 4 )*0. 

C YSM<5)*0. 

C YSMC 6 )*0 • 

YSMC 7 >*0. 

YSMC 8 )*0 • 

YSM9>*0. 

YSM( 10 )*0 * I 

C PLOT CONVOLUTED PROFILE 

CALL MQVEAC 0* »0, ) 

DO 940 M*1 9 10 

C DECONVOLUTE LASER FROM OBSERVED DATA ( /20 NORMALIZES ONLY) f 

940 YSMC M )*C YSMC M )/LC M ) >/20 5 

URITECSf*)' DECONVOLUTED COEFFS ARE' 

WRITE( 5f* )YSMC 1 ) ! 

URITEC 5 f # )YSM< 2) t 

URITEC 5 f # )YSM< 3) 

URITEC 5 f a ) YSMC 4) > 

URITEC 5f* )YSMC 5 ) 

URITEC 5f* )YSMC 6 ) 

URITEC 3 f # )YSM( 7) 

URITEC 5 f 4 )YSM< 8) 

URITEC 5f*)YSMC 9) 

URITEC 5 f a )YSM< 10 ) ! 

URITEC S f # )' THIS TEXT CAUSES LAST LINE TO PRINT 
C TO AVERAGE COEFFS WITH ANOTHER PROFILE 1112331121323 12211231131 

C Cl*8. 738/100 

C C2*4 . 383/100 l 

C C3*2. 346/100 

C YSMC 1 )*< YSMC 1 )+Cl )/2 

C YSMC 2 )*C YSMC 2 )+C2 )/2 

C YSMC 3 )*C YSMC 3 )+C3 )/2 

C URITEC 5f*)' AVERAGE COEFFS ARE ' 

C URITEC 5f# )YSMC 1 ) 

URITEC 5f* >YSMC 2 ) 

C URITEC 5*# )YSMC 3 ) 

C 311131132111111111111111111111111131111111211111111111111111311 

C ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦TO PLOT LASER PROFILE INSTEAD4#4 

C YSM( 1 >*.2433401 

C YSMC 2 >*.1676189 

C YSMC 3 )=. 1367337 
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ftu tk 


NT-; 


■*. V. 


I 


] 


c ♦*##***♦♦*****#*•«*♦**#**«#*♦**♦#*♦♦♦**♦♦*♦♦♦♦*♦ 

X-.01 

C FIRST PLOT IS NON NORMALIZED VERSION 

950 Y-YSMC 1 >4C0SC X/100 )4YSMC 2 >*COSC 24X/100 )+YSMC 3 >*CQSC 34X/100 > 
44YSMC 4 >4C0SC 44X/100 > 

4+YSMC 5 >*COSC 54X/100 ) 

4+YSMC 6 )4C0S< 64X/100 )+YSM< 7 )4CQSC 74X/100 HYSMC 8 >4CaSC 84X/100 > 
CC CALL DRAUAC X*Y> 

C URITEC 5* 4 )' X* Y' 

C WRITEC 5 * 4 )X* Y 

C FIND FIRST AND CENTRAL VALUES TO NORMALIZE IN STEP 980 
IF C X • NE » .01 ) GO TO 940 
FV*Y 

960 IF (X .NE. 315.01) GO TO 970 

C WRITEC 5*4 )'X IS 315 X AND Y ARE' 

C URITEC5*4)X*Y 

C WRITE< 5*4 >' 44*444444444444444' 

CV«Y 

C URITEC 5*4)' CENTRAL VALUE IS CV*' 

C WRITEC 5 * * )CV 

C WRITEC 5 *4 )' 44444*444*44444*44' 

970 X*X + 1 

IF ( X .LE. 629 ) GO TO 950 

URITEC 5*4 >' FIRST VALUE* FV* CENTRAL VALUE* CV' 

WRITEC 5*# )FV*CV 
C WRITEC 5*4 )' CV IS' 

C U)RITE(5 *4>CV 

WRITEC 5*4 )'THIS TEXT CAUSES LAST LINE TO PRINT' 

CO*C 15.6894164C < SRS/C 24PI ) >442 ) )/10000 
URITEC 5*4 >' CO' 

WRITEC 5*4 )C0 

URITEC 5*4 )' XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX' 


C 

CALL 

CALL 

MOVE ACC *0.) 
DRAWAC 629. *0. ) 

C 

CALL 

MOVEAC 0. * .2 ) 

C 

CALL 

DRAUAC 629. * .2 ) 

C 

CALL 

MOVEAC 0. * .4 ) 

c 

CALL 

DRAUAC 629* * .4 ) 

c 

CALL 

MOVEAC 0. * .5 ) 

c 

CALL 

DRAUAC 629. * .5 ) 

c 

CALL 

MOVEAC 0. * .6 ) 

c 

CALL 

DRAUAC 629. *.6) 

c 

CALL 

MOVEAC 0. * .8 ) 


C 300+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

C CALL DRAUAC 629 • * . 8 ) 

C CALL rt0VEA(0.*l. ) 

C CALL DRAUAC 629 • * 1 • ) 

C NOW PLOT NORMALIZED CONVOLUTED DATA 

CALL DU INDOC -3l5.*315.*0.*l. ) 

X*. 01-315. 

900 Y *C < YSMC 1 >4C0SC ( X/100 > HYSMC 2 )4C0SC 24X/100 ) 

4+YSMC 3 )4C0SC ( 34X/100 ) ) 
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*t-YSM< 4 >#COS< 4*X/100 ) 

*+YSM< 5 )*COS< 5*X/100 ) 

*4YSM< 6 )*COS< 6*X/100 ) 

*+YSM< 7 )#COS< < 7*X/ 100 ) >+YSM< 8 )*C0S< 8*X/100 ) 

*4YSM< 9 )#C0S< 9*X/lOO >+YSM< 10 )*CQS< 10*X/100 ) 

* )-CV >/( FV-CV ) 

CALL DRAUA(XfY) 

C WRITEC 5f*)X,Y 

X=X+1 

IF (X • LE ♦ 315) GO TO 900 
C PLOT THREE REFERENCE CURVES T=800f 1200t 1600 
CALL M0VEA<0 .f1* ) 

C Following is temperature of first reference curve 

T =800 

984 X=0 . 

CALL MOVEA< O.f 1. ) 

C Below is the theoretical gaussian fringe profile 

985 Y=EXP( -C0*< X*X )/T ) 

CALL DRAWA(XfY) 

X=X+10 

C In the following line 315 is used because pi radians of fringe plotted 

IF ( X .LE. 315 ) GO TO 985 
T =T +400 

IF (T .LE. 1600) GO TO 984 
C FIND SQUARE ERROR FOR T= 

C Start searching Tor correct temperature at the temperature below 

T=4rc 

WRI'E< 5f#)'TEMPf ERRORS 

986 £R=0 * 

X=0 . 

987 ER=£R+( i ( ( ( YSW 1 }*C0S< ( X/100 ) )+YSM( 2 >*C0S< 2#X/100 ) 

*FYSM< 3 )*C0S< < 3*X/100 > )FYSM< 4 )*C0S( 4*X/100 ) 

*+YSM< 5 )*C0S< ( 5*X/100 ) ;FYSM< 6 >*C0S< 6*X/100 ) 

#+YSM< 7 )*C0S< ( 7*X/ 100 ) )+YSM< 8 )*C0S< 8*X/100 ) 

*+YSM< 9 )*C0S< 9*X/100 )+YSM< 10 )*C0S< 10*X/100 ) 

* )-CV )/< FV-CV ) )- 

*< EXP( -C0*< X*X )/T ) ) )**2 ) 

X=X+10 

IF( X .LE* 315 ♦ ) GO TO 987 
C URITE(5f*)'TEMPr ERROR' 

URITE< 5f*)T,ER 

C WRITE< Sr#)' THIS TEXT CAUSES LAST LINE TO PRINT' 

T=T+10 

IF < T *LE* 1700 ) GO TO 986 
CL0SE<UNIT=3fDISP0SE='5AVE' ) 

C Next line applies to graphing on terminal 

CALL FINITT( 0f76C ) 

END 


.TYPE G. COM 

FORTRAN/LIST :FW2 : C . LST/SH0W:3 G. G 

linn g/linklibrary:fsrlid 

R G 
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. fortran/list:fu2:g.lst/shou:3 g.g 

* KR aN * 


. link g/linklibrary:fsplib 

TCxNK-W-Undef ined dlobalsl 
*V1RSZ 


• R G 

Zw ZS 

Httttttttt 

BACKGROUND INTENSITY IS 
980.7231 
75.53000 

NC BEFORE NORMALIZING AND AFTER SUBTRACTION OF BG 
41.25928 
136.7249 
11.27686 
120,3678 
72.98279 
119.8583 
125.0291 
413.9391 
433.7892 
816.1156 
1174.443 
1839.478 
2851.599 
4722.567 
6891.408 
8991.732 
10968.37 
11536.11 
11110.94 
9288.200 
6731.890 
4509.587 
2715.975 
1797.112 
918.6843 
663.2294 
319.3577 
169.7096 
303.6644 
149.3887 
141.9388 
114.2964 
104.1340 
-41.25934 
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UIX IS 

11536*11 

PRESSURE* COUNTS* Mi_ORAYLEIGHS* IN RAW DATA FORM 
LAST 2 COLUMNS ARE NEW PRESSURE AND NEW COUNTS 


PL 

CL 

KL 

NP NC 


957.1900 

584.0000 

0.5670000 

960.7752 

3.5765325E-03 

954*9500 

628.0000 

0.5580000 

958.5333 

1. 1851901E-02 

952.7000 

555.0000 

0.5550000 

956.2716 

9.7752654E-04 

950.3500 

610.0000 

0.5500000 

953.9160 

1.0434001E-02 

948.1900 

585.0000 

0.5510000 

951.7607 

6.3264635E-03 

946.200 0 

633.0000 

0.5710000 

949.7921 

1 . 0389834E-02 

943.850 0 

656.0000 

0.5890000 

947.4012 

1.0838059E-02 

941.6300 

843.0000 

0.6010000 

945.1664 

3.58B2030E-02 

939.2400 

872.0000 

0.6130000 

942.7611 

3.7602723E-02 

937.3100 

1119.000 

0.6200000 

940.8429 

7,07 44425E-02 

934.8000 

1339.000 

0.6190000 

938.3013 

0.1018058 

932.5200 

1745.000 

0.6170000 

936.0159 

0.1594539 

930.6400 

2381.000 

0.6200000 

934.1635 

0.2471890 

928.0900 

3541.000 

0.6200000 

931.5568 

0.4093725 

925.7500 

4799.000 

0.6090000 

929.2112 

0.5973770 

*23.3500 

6128.000 

0.6140000 

927.3290 

0.7794423 

921.5500 

7246.000 

0.6060000 

925.0110 

0.9507859 

919.3300 

7565.000 

0.6040000 

922.7765 

1.000000 

917.4700 

7163.000 

0.592000 0 

920.9374 

0.9631442 

914.9400 

6012.000 

0.5850000 

918.3515 

0.8051413 

912.7600 

4532.000 

0.5870000 

916.1726 

0.5835493 

910.6800 

3255.000 

0.5920000 

914.0910 

0.3909105 

*08.3000 

2208.000 

0.5960000 

911. 6835 

0.2354324 

906.3900 

1652.000 

0.5930000 

909.7880 

0 - 1557814 

904.0000 

1133.000 

0.5940000 

907.3646 

7 . 9635531E-02 

901.9500 

973.0000 

0.5890000 

905.3192 

5./491589E-02 

399.3400 

777.0000 

0.5940000 

902.6629 

2 • 7683303E-02 

397.0800 

680.0000 

0.5870000 

900.4006 

1 . 4711 161E-02 

395.0400 

747.0000 

0.5780000 

898.3924 

2.6322946E-02 

892.6400 

651 .0000 

0.5720000 

895.9683 

1 • 2949656E-02 

890.3300 

649.0000 

0.5740000 

893.6648 

1.2303872E-02 

888.2900 

621.0000 

0.5630000 

891 ,6171 

9.9077048E-03 

385.7800 

612.0000 

0.5600000 

889.0557 

9 ♦ 0267882E-03 

883.5200 

523.0000 

0.5520000 

886.7994 

-3.5765378E-03 

THIS TEXT CAUSES 

LAST LINE 

TO PRINT 




SRS 

75.53000 


YSM 

0.3696032 
0.2238235 
0. 1115526 
4 * 6207 257 E-02 
1.6699495E-02 
5.2267537E-03 
1 • 8393526E-03 
I.2564564E-03 
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ORIGINAL PACE IS 
OF POOR QUALITY 


3 . 5860834E-03 
3.0901073E-03 

BOU0PLOT OBSERVED DATA USING FOURIER COEFFS 


7 0 


DECONVOLUTED COEFFS ARE 
8.1958860E-O2 
6.9355413E-02 
5.2061040E-02 
2 . 4508 T 04E-02 
1.1238017E-02 
* . <i>828589E-03 
0 . OOOOOOO 
0.0000000 
0 . OOOOOOO 
0 . OOOOOOO 

THIS TEXT CAUSES LAST LINE TO PRINT 
FIRST VALUE. FV. CENTRAL VALUE. CV 
0.24380*9 -4.6711199E-02 

THIS TEXT CAUSES LAST LINE TO PRINT 
CO 

0.2267180 

ZXZBXXXRRRXKRaXKXBMRBRZBBXKXZXaXXBBXBEBEBBHXKIJLMwOPGxST > :*TaUWzXZC C \ I DF" 

‘GaHJbKLcNdOF’eRSfUVaWhYZiC 3a t >ZAKBClEFGImJ''niXJKMNOQRTUVXYZ\l 3t l&BKABDaEF lHIhJLaH 

OfFeQSdTcUWbXaY 'C\ ' 3"_ V 0A ; CCI" C ' HzEGaHxJwKLvNuDPtRsSTrUr*PPLQvSxTzU ! U X" aYdCaXJ 
tn_o40uBaC JD*aFeG#e+GjHo.JtNaL NUeOKPoRwS'TXeVlUlsYl Z&bC j 3st C d. AlBuC E(hFrGCI )eJp 
KnH >z,M*dNoOzG+eRpTCU»fVc*x: Y-hZs\ 3. jtv-G/bAn.BaDOdEpFC Hlf Ir J 3L2hMsN)P3hG3h-GrS}T4 
3 UpWzX5cY1Cu\" 36f _n . 8vA)C7dDKEaGwH ! I8bKf LKNoOrPvRxSC T JV^LI X8 .X9'ZC\8 T~_}/@CBxCu 
DrFoGK IfJaK7 ! MuNo0jGdR6 JSvUnVf W5"YuZlC c 34zT P_4p/_f 0A3 ! BrChE2 JFrHhl 1 JJoLf MOCNoPdQ/ 
xRroTaU. vVjX- Ys7h\» : 3«tf 10+zAoCdD+dlD*aEoGdH >zIoNeL( C MoOhP' "’QuSITcUtkCWrXJYbCXzXst 
l_e2BS"BwCaDKFeG* Ha.J#a2JLKoLjNe03P" TRaSuToVnWaXsZdC 3 3 • T ! _z3AxBvCtErFr 6AsBCtDuFG 
vHw.JKxM 1 ::6MaNzOQC R ! S 3UV"U Y " 'll a 3bTc_c7AdBeCEf FaHIhJLirtN aPQKRTUl VXYn.ZX 3tm80A" m8AB 
DEGHIKLIMOFQKST aUUX l YC h\ 3s_f 9BBeCHdFGcHbJK3LN ' 0 • PR W ST 3V ! UXCZCz J 1 =9 3 1 a_a : AxBCwEFG 
vIJuKMNOiORSUsVUYZ\3rtr;0ABDEFHp .BC9'.CCC8K/H6pV6p/V3u0C0oQ-ct*'lL' ijiYX12G9~T m ~3Ba 
O' 3::4.JuXsSiESr 60NCr7IVr8DQ_r9LZr JGU 1 p t Ur »BP 39 ' . L C 8r/H7 1 V5K0C2 ! QOt'T -n 1L#C Y< r2G&LT7. 
c3B*"0c3“o4JdX 1 35ExS ! x5Su60LNsC r 71 Vr 8BQ_r9LZr ! GUr » BP 39 ' . EC 8u/H7wV6t’0C4hU2 'T/tlL-i 
Y+a2G 
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.1 


...» . » fc-.-CTSCr 



i 



/ 

TEMP# ERROR 



400.0000 

0.2097060 


410.0000 

0.1945095 


420.0000 

0.1801474 


430.000 0 

0.1665846 


440.0000 

0.1537879 


450.0000 

0.1417260 


460.0000 

0.1303690 


470.0000 

0.1196891 

/ 

480.0000 

0.1096597 


490.0000 

0. 1002555 

; 

500.0000 

9. 1452442E-02 


510.0000 

8. 3228029E-02 


520.0000 

7.5560570E -02 


530.0000 

6.8429396E-02 


J40.0000 

6. 1815027E-02 


550.0000 

5.5698719E-02 

- 

560.0000 

5.0062709E-02 


570.0000 

4. 4890013E-02 


580.0000 

4 . 0164385E-02 

t 

590.0000 

3.5870302E-02 


600 • 0000 

3. 1992938E-02 

- 

610.0000 

2.8518166E-02 


o20 ♦ 0000 

2 . 5432343E-02 


630.0000 

2.2722503E-02 

* 

640.0000 

2. 0376189E-02 

- 

650.0000 

1.8381432E-02 

ji 

660 • 0000 

1 . 6726822E-02 


670.0000 

1.S401303E-02 

' 

680.0000 

1 .4C;94365E— ,)2 


690.0000 

1.369585 IE -02 


700.0000 

1.3296O20E-O2 


710.0000 

1 . 3 18550oE-v>2 


720.0000 

1 . 3355302E-02 


730.0000 

1 . 37O6733E-02 


740.0000 

1.4501 466E-02 


750.0000 

1.546i.o8E-02 


760.0000 

1 . 6668990E-02 


770.0000 

1.811&590E-02 


780.0000 

1 .9797074E-02 


790.0000 

2. 1703515E-02 


800.0000 

2 . 3829229E-02 

- 

810.0000 

2.6167762E-02 


320.0000 

2.8712889E-02 


830.0000 

3. 1458605E-02 


840.0000 

3.4399092E-02 


850.0000 

3.7528761E-02 


860.0000 

4.0842198E-02 


870 . 0000 

4.4334162E-02 


TC880.0000 

4.7999568E-02 


\ 

I 

i 
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TYPE 0*0 


C PROG TO READ AND MODIFY LASER PROFILES 
C TEMP< LASER) PRESSURE( LASER ) COUNTS< LASER ) 

DIMENSION YSM< 35 ) t YSCM< 35 >• YSSM< 35) 

REAL TL( 130 ) 

REAL PL( 130 ) 

REAL CL< 130 ) 

REAL NP< 130 ) 

REAL NC( 68 ) 

OPEN< UNI T=3 f NAME* ' FU2 • FP1 1 * DAT' r READONLY f TYPE®' OLD' ) 
WRITE<5 f*)' TL NP CL 

DO 100 1*1 , 68 

READC 3? 50 ) TL( I )* PL( I )fCL( I ) 

IF ( I «LE* 4 ) GO TO 100 
C START TEMPERATURE COMPENSATION 

C NEW PRESSURE < NP )=PRESSURE OF LASER< PL >*TEMP FACTOR 

NP( I )=PL( I )*< 298 ♦ /< 273 * f TL( I ) ) ) 

C NORMALIZE THE COUNTS 

NC< I )=CL( I )/32350 

URITE( 5t * )TL< I ) t NP< I )»CL( I >*NC( I ) 

50 FORMAT ( 9X, F6. 3, F7 . 2t 15X, F6 . 0 ) 

100 CONTINUE 

C NOW SMOOTH THE LASER SCAN USING FQURIER/RQBLE METHOD 
PI=3 * 1415926535 
°ITE( 5* # >' NP AT 1 * 2 ? 3*4*5' 
ulRITE( 5?* )N P( 1 ) 

WRITE( 5*# )NP< 2 ) 

WRITE( 5» * )NP< 3 ) 

URITE<5**)NP< 4) 

URITEC 5* 4 )NP< 5 ) 

WRITE( 5*# )'NP AT 68t69' 

WRITEC 5* * )NP< 68 ) 

WRITE( 5*# )NP< 69) 

NP< 1 )=NP( 5 ) 

NP( 20 )=NP( 68 ) 

NP( 2 )-NP( 17 ) 

NP( 3 )=NP( 19 ) 

NP( 4 )=NP( 21 ) 

NP< 5 )*NP< 23 ) 

NP( 6 )=NP( 25 ) 

NP( 7 )=NP( 27 ) 

NP< 8 )=NP( 29 ) 

NP< 9 )=NP< 31 ) 

NP( 10 )=NP( 33 ) 

NP< 1 1 ")=NP( 35) 

NP< 12 )=NP( 37 ) 
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NPC 13 )=NPC 39) 

NPC 14 )*NP( 41 ) 

NPC 15 )=NPC 43 ) 

NP( 16 )*NPC 45 ) 

NPC 17 )*NPC 47) 

NPC 18 )-NPC 49 ) 

NPC 19 >=NP< 51 ) 

W*1 

WRITEC 5t * )' NEW PRESSURE VALUES HP' 

103 WRITEC 5f*)WfNPCU) 

W*W41 • 

IF CU .LE. 20) GO TO 103 
DO 106 M»1 »21 

C FSR IN PRESSURE UNITS 

C SRS=53.2464 

C SRS«43 • 3838 

SRS*75.5332 
SUNC s 0 • 

SUNS=0 * 

C2=6.28318*M/SRS 
DO 107 I-5»68 
IF ( I .NE. 5 ) GO TO 104 
OPS-NPC 141 )-NPC I ) 

GO TO 105 

104 CONTINUE 
DPS-NPC I >-NP< I-t > 

105 CONTINUE 

IF ( M .NE* 1 ) GO TO 116 
CC WRITEC 5f* )DPS 

CC WRITE' 5f* )NP< I ) 

116 CONTINUE 

SUNOSUNCK NCC I >/M )*C 2.*SIN< C2*NP< I ) >#< C SIN< C2*0PS/2. > '**2 

* >+COS< C2*NPC I ) )*SINC C2#DPS ) ) 

SUNS=SUNS4( -NCC I )/M >*< 2.*C0S< C2*NPC I ) )*< < SIN< C2*DPS/2. ) >**2 

* )-SIN< C2 *NPC I ) )*SINC C2KDPS ) ) 

107 CONTINUE 

YSMC M >=SGRT< SUNS*SUNS4SUNC*SUNC )/PI 
YSCMC M )-SUNC/PI 
Y5SHC M )=SUNS/PI 

106 CONTINUE 
WRITEC 5 t*Y YSM' 

WRITEC 5? # )YSM( 1 ) 

WRITE( 5f* )YSM< 2 ) 

WRITEC 5f* )YSMC 3 ) 

WRITEC 5f* )YSMC 4 ) 

WRITEC 5f* )YSHC 5) 

WRITEC 5f# )YSMC 6) 

WRITEC )YSM( 7 ) 

WRITEC 5f* )YSMC 8 ) 

WRITEC 5f* )YSMC 9) 


5 *^ 
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u-'r. 




WRITE* 5t* )YSM< 10 ) 

WRITE* 5 1 * )YSh( 11 ) 

UIRITE( 5f* )YSM( 12 ) 

WRITE< 5t* )YSM* 13 ) 

WRITE* 5t* >YSM< 14 ) 

WRITE* 5t*)YSM* 15) 

WRITE* 5t* >YSM< 16 ) 

WRITE* 5t* )YSM< 17 ) 

U)RITE( St* >YSM< 18 ) 

WRITE* 5 1 # >YSM< 19 ) 

URITE( St* >YSM* 20 ) 

WRITE( St* )YSh< 21 ) 

CALL INITTC30) 

CALL DWINDQ<0.»62?*t-l.fl. ) 

CALL TWINDQ* SO 1 900 1 50 1 800 ) 

CALL rtOOEA* 0. tO. ) 

X=.01 

900 Y=YSM( 1 )*COS< X/100 >fYSM* 2 )*COS< 2*X/100 >+YSM* 3 )*COS( 3*X/100 )+YSM( 4 ) 
**COS< 4*X/ 100 )+YSM< 5 >*COS< 5*X/100 >+YSM* 6 >*COS< 6*X/100 HYSH* 7 )*COS 
*( 7*X/ 100 > 

*4YSM( 8 )*COS( 8*X/100 ) 

* + YSM( 9 >*CQS< 9*X/ 100 )+YSM( 10 )*COS 
*< 10*X/100 ) 

CC *+YSM* 11 >*COS( 11*X/100 HYSM* 12 >*COS( 12*X/100 ) 

C *+YSM* 13 )*COS( 13*X/100 )+YSM( 14 )*COS( 14*X/100 )+YSh( 15 )*COSC 15*X/100 ) 

C *4YSM* 16 )*CQS< 16*X/1C0 HYSM* 17 )*COS* 17*X/100 )+YSM( 18 >*CC3< 18*X/100 ) 

C *+YSM< 19 )*COS( 19*X/ 100 HYSM* 20 )*COS( 20*X/100 >4YSM* 21 >*COS( 21*X/100 ) 

C IF ( X . GE. 2 ) GO TO 910 

WRITE* St* >'Xt Y' 

WRITE* 5* * >X t Y 
910 CALL DRAUA*XrY> 

X=X+1 

IF (X .LE. 629 ) GO TO 900 
CLOSE* UNIT=3t DISPOSE-' SAOE' ) 

CALL FINITT* 0 1 760 ) 

END 


♦TYPE 0 ♦ COM 

FORTRAN/LIST ; FW2 : 0 . LST/SHOU 13 0.0 

LINK o/linklibrary:fsplib 

R 0 


.60 
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* FORTRAN/LIST : FW2J 0. LST/SHOU 53 0.0 

.MAIN. 

. LINK o/linklibrary:fsplib 

?LINK-W-Undef ined globelsi 
SOIRSZ 

. R 0 


TL NP 



CL 

NC 

?Err 5 

in routine “.MAIN." 

line 

10 



?Err 5 

in routine " .MAIN.' 1 

line 

10 




23.77200 

939.0898 

3509.000 

0.1084*99 

23*88800 

940.5798 

3409.000 

0.1053 787 

23*99800 

941.9172 

3434.000 

0.1061515 

24*10500 

943.1827 

3317.000 

0.1025348 

24*21600 

944*3544 

3374.000 

0.1042968 

24*31800 

945.4337 

3304.000 

0.1021329 

24*44800 

946.3330 

3333.000 

0.1030294 

24.64100 

946.7205 

3189.000 

9 . 8578051E-0 

24*72700 

947.6982 

3211.000 

9 . 92581 17E-0! 

24*77300 

948.9128 

3386.000 

0.10466 77 

24*81400 

950.1931 

3363 ♦ 000 

0.1039567 

24.84900 

951*4221 

3287.000 

0. 1016074 

24*88700 

952.6312 

3385.000 

0.1046368 

24*97300 

953.2164 

3454.000 

0.1067697 

24*99300 

954.5724 

3323.000 

0.1027202 

25*02300 

955.7662 

3642.000 

0.1125811 

25.05100 

956.8363 

3622.000 

0.1119629 

25*06800 

957.8914 

3722.000 

0.1150541 

25*07500 

959.3185 

3837.000 

0. 1186090 

25*08100 

960.7289 

4064.000 

0.1256260 

25.10700 

961.6247 

4104.000 

0. I268c24 

25.11100 

963.1113 

4476.000 

0.1383617 

25.11600 

964.4546 

4941.000 

0.1527357 

25.12700 

965.7284 

5274.000 

0.1630294 

25.12900 

967.0714 

5498.000 

0.1699536 

25.13400 

968.3645 

6284.000 

0.1942504 

25. 14800 

969.1286 

6991.000 

0.2161051 

25.14300 

970.1844 

8488.000 

0.2623802 

25.13200 

971.5097 

11030.00 

0.3409583 

25.12700 

972.7354 

16357.00 

0.5056260 

25.12600 

974.1201 

24998.00 

0.7727357 

25.12400 

975.5840 

32280.00 

0.9978362 

25.12700 

976.8936 

32312.00 

0.9988254 
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25 * 14000 

9 77.7606 

28974.00 

0.8956414 

25.14000 

979.2099 

22393.00 

0.6922102 

25.14600 

980.5896 

16566.00 

0.5120866 

25.15100 

981.8824 

12226.00 

0.3779289 

25.15600 

983.2053 

9171.000 

0.2834930 

25.15900 

984.4948 

7474.000 

0.2310355 

25.16800 

985.0446 

6617.000 

0.2045441 

25.16500 

986.5237 

5754.000 

0.1778671 

25.16200 

987.8430 

5349.000 

0.1653478 

25.15800 

988.9757 

4901.000 

0.1514992 

25.15800 

990.1851 

4684.000 

0.1447913 

25.16000 

991.4576 

4362.000 

0.1348377 

25.16000 

992.7970 

4224.000 

0.1305719 

25.15200 

993.5732 

4035.000 

0* 1247295 

25.15100 

994.8059 

3970.000 

0.1227202 

25.14900 

996.0220 

3882.000 

0.1200000 

25.14500 

997.3847 

3711.000 

0.1147141 

25.14300 

998.7208 

3633.000 

0.1123029 

25.14100 

999.9669 

3773.000 

0.1166306 

25.14200 

1000.813 

3606.000 

0. 1114683 

25.13800 

1002.316 

3611.000 

0.1116229 

25.13900 

1003.622 

3613.000 

0. 1116847 

25.14200 

1004.941 

3541,000 

0.1094590 

25.14500 

1006.290 

3732.000 

0.1153632 

25.15300 

1007.503 

3567 . 000 

0.1102628 

25.15700 

1008.749 

3707.000 

0,1145904 

25.15800 

1 009 . 565 

3607,000 

0.1114992 

25.15800 

1011.004 

3609.000 

0.1115611 

25.16200 

1012.160 

3605.000 

0.1114374 

25.15600 

1013.529 

3551.000 

0.1097682 

25.14300 

1014.623 

3434.000 

0.1061515 

NP AT 1 r 2 f 3 1 4 

f 5 



0.0000000 




0.0000000 




0.0000000 




0.0000000 




939.0898 




NP AT 68*69 




1014.623 




0.0000000 




NEW PRESSURE 

VALUES NP 



1.000000 

939.0898 



2.000000 

952.6312 



3.000000 

954.5724 



4.000000 

956.8363 



5.000000 

959.3185 



6.000000 

961.6247 



7.000000 

964.4546 
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8.000000 

967.0714 


9.000000 

969.1286 


10.00000 

971.5097 


11.00000 

974.1281 


* 

1?. 00000 

976.8936 


13.00000 

979.2099 


14.00000 

981.8824 


13.00000 

984.4948 


16.00000 

986.5237 


17.00000 

988.9757 


18.00000 

991.4576 


19.00000 

993.5732 


20.00000 

1014.623 


YSM 



0.2254809 



0.1613598 



0.1071364 



9.4267033E-02 



7 • 42991 19E-02 



5.5807292E-02 



c 

2. 96831 63E -02 



2.8193893E-02 



1.8680153E-02 


j 

1 • 8559434E-02 


1 

4 • 4040936E-03 



9. 1322120E-03 



4.7645955E-03 



7.8142360E-03 



4.0516402E-03 


: 

2 • 6422290E-03 



5.1584304E-04 



2 . 9768483E- 03 



3 • 2191928E-03 




3.7323066E-03 




1 .9601721E-03 




XfY 




9.9999998E-03 

0.8134672 


Xf Y 



t .010000 

0.8127617 


Xf Y 



! 

2.010000 

0.8106764 


X 

Xf Y 



X 

3.010000 

0.8072207 

i 

» 

1 

Xf Y 


! ’v ' 

4*010000 

0.8024113 

- 


Xf Y 



5.010000 

0.7962701 

1 

XfY 


j 

6.010000 

0.7888262 

* t 

XfY 



< 1 

7.010000 

0.7801136 


XfY 


1 

8.010000 

0.7701731 


Xf Y 


1 

9.010000 

0. 7590501 


XfY 


i 

10.01000 

0.7467955 


1 I 

XfY 


; I 

11.01000 

0.7 334650 


1 

( i 

XfY 



1 

7 te. oiooo 

0.7191188 


f 1 


I 
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APPENDIX E, SHADOW HEIGHT 



Shadow Height in Two Dimensions 

To explain the concept of shadow height, this section shows how shadow 
height is calculated for the simplified case of two dimensions. 

The first step is to find: 

M » distance from observer to shadow along line of observations from 
a and 0, where 

a « angle from local vertical to line of observation, where positive 
values are toward the Sun, and 

0 ** change in latitude between observer and sunset as shown. 

M is found in terms of a and 0. Then, shadow height will be found in 

terms of M and a. 

Let 0 = a - 0 , h = M cos 0 

h ^ h m r - r cosQ 

cos 0 cos (a-0) co f (a-0) 

which applies for a±90°. 

Shadow height (SH) is determined from M and a (see Figure E-l), using the 
cosine law, where Y = 90 - |a| 

(SH +i) 2 ®r 2 +M 2 - 2rM cos (90 + Y) 

SH + r = [r 2 + M 2 - 2rM cos (90 + 

SH ■■ {r 2 + M 2 - 2rM cos (90 + Y]* - r 

SH » (r 2 + M 2 - 2rM cos |a|]^ - r 
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because cos (90 + ¥) « cos (90 + (90 - |a|)) 


■ cos (180 - |a|) 

■ -cos |a| 

Shadow Height Program Rotations Used 

Two shadow height programs are included. One of these programs (L.SH) 
is a simple version, but it requires the solar depression angle Tl, and the 
azmuth (or bearing, i.e., E from N) angle from the Sun of the observation 
be input (in degrees). These two angles may not be known, but they are not 
required in the more complex version of the program (M.SH). Instead, the 
time and place are Input and these two angles are calculated in M.SH. The 
rest of M.SH is the same as program L.SH. 

The components of the vector from the F.arth to the Sun in "inertial 
space" are found in M.SH. This vector is called VI. By performing a 
series of rotations, this vector is transformed into the frame of the 
observer, so the resulting vector VL gives the direction to the Sun from 
the observer in terms of local zenith and azmuth. 

To demonstrate this, an example is shown below. In Figure E-2 the 
components of VI are X, Y, and Z. With an angle of 66° for these 
components are 0, 0.9, and 0.4, respectively. The next rotation produces 
the "Geographic Vector," VG, which is fixed with respect to Earth and 
points toward the Sun. The amount of rotation, -GSTR, depends on the 
time. 

The next rotation is the longitude rotation. If the observer is at 
the 0° longitude, in England, as he is in this example, the rotation Is 0. 

The next rotation is the latitude rotation. Here, the coordinate 
system is rotated 90° about Y; then the coordinate system is rotated 
through an angle that equals latitude, so the Z axis is in the local 
vertical. The X axis points south and the Y axis points east. As in the 
other rotations, the new components for the unit vector pointing to the Sun 
In this rotated coordinate frame are computed. These new vector components 
are VL(1), VL(2), and VL(3). 

Finally, the observer has the components of a unit vector pointing 
toward the Sun in terms of local zenith, south, and east, so he can easily 
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calculate the zenith and bearing angles to the Sun. The zenith angle; 
l.e., the angle from the zenith to the Sun is 


■ tan 


-1 V y 2 + x^~ 


which can be visualized by translating the XY plane component along Z to 
the maximum Z value. 

To find the bearing, or azimuth angle, 


6 “ tan 


-1 


1 

x 


is computed, where 6 is measured to the E from S. To convert this to a 
"bearing" angle which is defined as angle E of N, 


Azimuth * Bearing ■ 180 - 6 


TABLE E-l. SHADOW HEIGHT PROGRAM VECTOR COMPONENTS 



X 

Y 

Z 


VI 

0 

0.9 

0.4 

inertial vector 

VG 

0.9 

0 

0.4 

geographic vector 

VJ 

0.9 

0 

0.4 

after longitude rotation 

VL 

0.05 

0 

0.95 

after latitude rotation 
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Shadow Height Program s 


TYPE M.SH 

C ni.sh SHADOW HEIGHT PROG. BASED ON LOCAL TIME 

DIMENSION VI< 3 >fB V< 3 >fVG< 3 )fVLC 3 >f VJ< 3 >fPV< 3 >f ZZ< 3f3 ) 

DIMENSION C(3f3)fCT< 3f3> 

CALL INITT< 30 > 

CALL DUIINDO( 7200 f67000. fO. f250. ) 

CALL TWIND0<50f900f50fQ00) 

CALL MOVEA< 0. fO. ) 

CALL MOVE A( 3600. fO. ) 

CALL DRAUA< 3600 . f 10 . ) 

CALL MOVEA( 7200. fO. ) 

CALL DRAWA( 7200. fO* ) 

RAD=57. 295779 

C LATITUDE AND LONGITUDE OF OBSERVER IS 

ALA=64. 86001/RAD 
ALQ--147 .84711/RAD 
C ALA=49./RAD 
C AL0=-122. /RAD 
1 DAY -61 
C I DAY =46 
I YR= 1984 

3 SECS-O » 

CALL M0VEA(0.f0.) 

C 4 WRITE( 5f#)' ENTER YEAR < INTEGER )f DAY ( INTEGER )f SECONDS < REAL >' 

4 MS=1 

C #**##**#«#*#**** 

c tttttttttttttrt 

c 5 READ<5**) IYRfIDAYfSECS 

IF ( IYR.LT. 1901. OR. IYR.GT. 2099 ) STOP 

CALL SUN( IYRfIDAYfSECSfGSTfSLONGfSRASNfSDEC > 

C WRITEC 5 f 200 ) 

C WRITE( 5 f 300 ) IYRf IDAY » SECSf GST ? SLONGf SDEC 

C GO TO 5 

COO FORMATf 214 f FI 0.2) 

200 FORMAT( 8Xf ' IYR' f 6X f ' IDAY' f6Xf ' SECS' f7Xf ' GST' f5Xf' SLONG' f 

#5Xf ' SDEC' ) 

300 FORMAT( IX f 2I10f4F10 . 3 ) 

C *************JM********************************************* 

C NOW FIND INERTIAL COORDS (VI) FROM SDEC AND SLONG 

P»< ( 3.141592/2 >-< SDEC/RAD ) ) 

SLONG=SLONG/RAD 
VI1-SIN( P >#CQS< SLONG ) 

VI2=< SIN< P >*SIN< SLONG ) ) 

VI3=C0S(P> 

VI( 1 )=VI1 
VI( 2 >=VI2 
VI( 3 )=VI3 

C WRITE< 5 f* )' INERTIAL VECTOR VI X Y Z COMPS ARE' 

C WRITE(5f*)VI 

C NOW TRANSFORM THE INERTAIL COORDINATES (VI) TO GEOGRAPHIC ( VG ) 

GSTR»GST/RAD 
CALL ROTXYZ( -GSTRfZZf3) 
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CALL TRF< VI rZZ» VG > 

C URITEC 5t*>' GEOGRAPHIC VECTOR VG X Y Z COMPS TO SUN FROM EARTH' 

C WRITEC 5t # )VG 

C «»««********«**«]K]|C*«*«******)K*****¥: ******«****»«***»««******«]((***«« 

C NOW TRANSFORM THE GEOGRAPHIC COORDS < VG ) TO LOCAL AZIMUTH AND ZENITH 

C NEXT ROTATE LONG ABOUT Z AXIS 
XI=90./RAD 

CALL ROTXYZC -AL0tZZ*3) 

CALL TRFC VG*ZZ>VJ > 

C URITEC 5r* )' VJ IS' 

C URITEC 5>*)V.J 

C NEXT ROTATE LAT t AND 90 DEGREES ABOUT Y 
E=-ALA+XI 

CALL ROTXYZC E»ZZ* 2 ) 

CALL TRFC VJ>ZZ,VL ) 

C URITEC 5*# )' ZENITH FROM ACOS Z TO SUN IS ' 

X=VL< 1 > 

Y=VL( 2 ) 

Z=VL< 3 ) 

SQ=SQRTC C X**2 )+< Y**2 ) ) 

AZ=< ATAN2C 5Q>Z ) )*RAD 
C URITEC 5t*)AZ 

TH=( ATAN2C YfX ) >*RAD 
B=180-TH 

C URITEC 5*# )' AZIMUTH TH IS' 

C URITEC 5f*)TH 

C URITEC )' BEARING TO SUN IS' 

C WRITEC 5* # )B 

C **# NOU FIND SHADOW HEIGH! ***«#M*«********************** 

C FIND SOLAR DEPRESSION ANGLE 

T I=AZ-°0 » 

C TI = 45 » 

C WRITEC 5t* V SOLAR DEPRESSION ANGLE IS' 

C WRITEC 5 f # )TI 

C INPUT GEOGRAPHIC BEARING ANGLE OF OSERVATION 

C URITEC 5r* >' INPUT GEOGRAPHIC BEARING ANGLE OF DBS (DEG)' 

C REA DC 5 r # )GB 

GB=*01 

COO URITEC 5f* )' GEOGRAPHIC BEARING ANGLE OF OBS IS' 

C URITEC 5t* )GB 

900 MS2=1 

C INPUT ENITH OF OBSERVATION ( DEGREES) 

C URITEC 5?# )' INPUT ZENITH OF OBS' 

C READ(5f*)ZI 

ZI=60 ♦ 

C WRITEC 5r* )' ZENITH OF OBSERVATION IS ' 

C URITEC 5 f * )ZI 

C CALCULATE BEARING ANGLE OF OBSERVATION FROM SUN 

910 SI=GB-B 

C SI = 45 • 

C URITEC 5* * )' BEARING ANGLE OF OBSERVATION FROM SUN IS' 

C URITEC 5>* )?T 


53 



ooooo ooo 





NOW RUN SHADOW HEIGHT PROG L.SH 

Z*ZENITH ANGLE* S=»AZIMUTH, T*SOLAR DEPRESSION ANGLE 
RADIUS OF EARTH IS KM IS: 

RA*6378. 

THE FOLLOWING 3 VALUES ARE INPUTS TO PROG (ABOVE) 

ZI IS ZENITH INPUT (DEGREES) OF OBSERVATION 

SI IS AZIMUTH INPUT (DEGREES) BEARING ANGLE FROM SUN OF OBSERVATION 
TI IS SOLAR DEPRESSION ANGLE INPUT (DEGREES) 

NOW CONVERT THESE TO RADIANS: 

Z=ZI/RAD 
S=SI/RAD 
T=TI/RAD 

CI = l/( ( SIN( T ) >*#2 ) 

C2=( ( ( COS( Z ) )**2 )#( ( COS( T ) )*#2 ))/((( SIN( Z ) )*#2 >*( ( SIN( T ) )#*4 ) ) 

C3=( 2#C0S( Z )*COS( T )*COS( S ) )/( SIN( Z )#( ( SIN( T ) )**3 ) ) 

C4=( COS( T >*COS( S ) )/SIN( T ) 

C5=C0S( Z )/( ( ( SIN( T ) )##2 )*SIN( Z ) ) 

C6=C4+C5+( ( C1+C2+C3 )## . 5 ) , 

AP=( SIN( Z )#C6 )-( COS( Z ) ) ' 

p= RA/AP 

AG=( ( ( RA*#2 )+( P**2 )+( 2*RA*P*( COS( Z )>>)*#. 5 )-RA 
C URITE( 5,#)' HEIGHT ABOVE GROUND IS' 

C WRITE(5f#)AG 

GB=GB+90. 

c write: 5»# )' =====================' 

IF (GB .EO. 90.01) GO TO 950 

IF (GB .EO. 270.01) GO TO 955 i 

IF (GB .EO. 360.01) GO TO 960 i 

IF ( GB .EO. 450.01) GO TO 97 0 j 

GO TO 1000 

950 32=AG ’ 

GO TO 1010 : 

955 31=AG 

GO TO 1010 
960 W1=AG 

GO TO 1010 
970 S4=AG 

GO TO 1020 ; 

1000 S3=AG .V 

1010 IF( GB .LT. 280 ) GO TO 900 

C NOW ZENITH OBSERVATION 

ZI=0.01 
GO TO 910 

C WRITE IN ORDER OF N E S W ZENITH i 

C020 WRITE( 5f#)S2»S3»Sl»Ul*S4 

1020 LL=1 . 

IF ( SECS .GT. 61200. ) GO fO 2000 
IF (SECS .LT. 10800.) GO TO 2000 
CALL DRAWA (SECS»S4) 

2000 SECS=SECS+1200. 

LF( SECS .LT. 36400. i GO TO 4 

IDAY = IDAV +8 ! 

C URITE( )' DAY' 


U 
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c 


.TYPE 


.TYPE 


URITE< 5#* )IDAY 
GO TO 3 

CALL FINITTC 0,760 > 
END 


H.COM 

fortran/list :FW 2:rt.LST/SH0w : 3 n.sh 
link rtf allii/libraky:fsplib 
r rt 


SUN. CAL 

SUBROUTINE SUN( I YR f IDAY , SECSf GST f SLONG* SRASN * SDEC ) 

DATA RAD /57. 29578/ 

DOUBLE PRECISION DJ»FDAY 

IF( IYR.LT. 1901. OR. IYR.GT. 2099 ) RETURN 

FDAY=SECS/96400. 

0J = 365*( IYR-1900 )+( IYR-1901 )/4+IDAY+FDAY-0 .5D0 
T=D*J/36525. 

UL=DMOD< 279. 696678+0. ?8564^3354#DJ, 360. DO ) 

GST~DrtOD< 279.c»90983+0.9856473354*DJ+360.*FDAY + l80. , 360. DO ) 

G=DMUDC 358.475845+0. 985600267*DJ* 360. DO)/RAD 

SLONG-UL+< 1. 9194^-0. 004789*T )*EIN( G >+0 . 020094*SIN< 2.*G ) 

l)BLIQ~< 23. 45229-0. 0130125*T )/RAD 

3LP=< SLONG-O. 005686 VRAD 

SIND=SIN( OBLIQ )*SIN< SLP ) 

COSD=SQRT( 1 • ~*SIND#*2 ) 

SDEC=RAD*ATAN< 3IND/C0SD > 

COT AN=COS( OBLIQ >/SlN( OBLlQ ) 

3RASN=t80 • ~RAD#ATAN2< COTAN*SlNO/COSDf -COSC SLP )/COSD ) 

RETURN 

END 
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SUBROUTINE TRF<X*A*XT) 
DIMENSION X< 3 ) * A( 3 * 3 ) * XT< 3) 
DO 1 1 = 1*3 
XT( I )=0.0 
DO 1 J=l*3 

XT( I )=XT( I >+A< I * J >*X< J ) 

1 CONTINUE 

RETURN 
END 

SUBROUTINE ROTXYZ< A*B* IROT > 
DIMENSION B< 3 * 3 ) 

A1=CQS( A) 

A2=SIN< A) 

GO TO ( 1*2*3)* IROT 
C ROTATION ABOUT THE X-AXIS 

1 B( 1*1 >=1.0 

B( 1*2)=0.0 
B< 1*3 )=0 .0 
B< 2*1 )=0 • 0 
B( 2*2 >=A1 
B( 2 * 3 )=-A2 
B( 3 * 1 >=0 » 0 
B< 3*2 )=A2 
B< 3*3 >=A1 
RETURN 

q Bp9If^8¥ A&0UT THE Y-axis 

B( 1*2 >=0 ♦ 0 
B ( 1*3 )=-A2 
B( 2*1 >=0.0 
B( 2 * 2 ) =1 A 
B( 2*3 )=0.0 
B< 3*1 )=A2 
B( 3*2 )=0 . 0 
B( 3*3 >=A1 
RETURN 

ROTATION ABOUT THE Z-AXIS 
B( 1*1 )=A1 
B( 1*2 >=-A2 
B< 1*3 >=0.0 
B( 2*1 >=A2 
B( 2*2 >=A1 
B< 2*3 >=0.0 
B( 3*1 >=0.0 
B< 3*2 )=0 *0 
B( 3*3 )=1 .0 
RETURN 
END 
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ORIGINAL PAGE IS 
OP POOR QUALITY 

APPENDIX F. EXAMPLE OF RAW DATA 


Scan 





Encoder 


Counts 

Photometer 




Mirror 

Number 

T,° 

C Pressure 

Position 






Position 

1 

14 

24 

246 

697. 

76 

3876 

0 

456 

0 206 

0 

0 

0000 0 0 

1 

3 

35 

24 

248 

699 

86 

3969 

-3 

476 

0 206 

0 

0. 

0000 0 0 

1 

3 

36 

24 

245 

902 

33 

4056 

0 

530 

0 203 

0 

0 

0000 0 0 

1 

3 

3 7 

24 

246 

904 

17 

4150 

143 

691 

0 211 

0 

0 

0000 0 0 

1 

3 

3 B 

24 

251 

906 

65 

4240 

0 

763 

0 203 

0 

0 

0000 0 0 

1 

3 

39 

24 

240 

909 

46 

4330 

9 

1067 

0 205 

0 

0 

0000 0 0 

1 

2 

4) 

24 

247 

911 

79 

4420 

6 

1347 

0 210 

0 

0 

0000 0 0 

1 

3 

41 

24 

254 

914 

02 

4510 

3 

1664 

0 213 

0. 

0 

0000 0 0 

1 

3 

42 

24 

251 

916 

51 

4600 

-3 

2013 

0 214 

0 

0 

0000 0 0 

1 

3 

43 

24 

250 

918 

82 

4690 

3 

2172 

0 212 

0 

0 

0000 0 0 

1 

3 

14 

24 

250 

920 

69 

4778 

-3 

2065 

0 211 

0 

0 

0000 0 0 

1 

3 

45 

24 

240 

923 

32 

4664 

-3 

1754 

0 215 

0 

0 

0000 0 0 

\ 

•» 

46 

24 

2 ? 9 

925 

82 

4932 

12 

1438 

0 215 

0 

0 

OOCO 0 0 

1 

3 

<7 

24 

242 

628 

33 

5040 

137 

1096 

0 216 

0 

0 

0000 0 0 

1 

3 

46 

24 

250 

930 

42 

5124 

0 

835 

0 213 

0 

0 

0000 0 0 

1 

3 

49 

24 

243 

932 

98 

5210 

6 

656 

0 214 

0 

0 

0000 0 0 

1 


50 

24 

246 

935 

17 

5292 

6 

62? 

0 212 

0 

0 

0000 0 0 

1 

3 

5i 

24 

24 1 

937 

le 

5376 

-3 

496 

0 211 

0 

0 

0000 0 0 

1 

3 

52 

24 

242 

939 

55 

5458 

-3 

514 

0 215 

0 

0 

0000 0 0 

1 

3 

53 

24 

235 

942 

09 

5542 

3 

462 

0 211 

0 

0 

0000 0 0 

; 

3 

5 4 

24 

240 

944 

49 

5624 

-3 

437 

0 210 

0 

0 

0000 0 0 

1 

3 

55 

24 

239 

946 

89 

5704 

15 
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APPENDIX G. VERIFICATION OF EQUATION (A-l) , APPENDIX A 


We can find the temperature T In terms of width by substituting 1 I 
for I In equation (A-l) 


max 


jM _ 2 ) h 


where X max = \ 2«kTA 


Substituting Into equation (A-l) gives: 




h (M c 2 \H 


i 2kl 


(AA) ‘ 


exp >2m*“ 


Solving for T, we obtain: 


In Os) = 


-2.196825683 x IQ 26 


T = 3.169349519 x 10 26 (AA) 2 


where AX is the HWHM in meters, 
o 

To input A rather than m, we have: 

T = 3.169349519 x 10 6 (AA) 2 
o 

where AX is the HWHM specified in A. 

To put the formula in terms of FWHM, multiply the numerical constant by 

(i) 2 * i 
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T = 7.92337 x 10 5 (AX) 2 

, O 

where AA is the FWHM in A. 

This is consistent with the result of Wark [3] who obtained T « 7.89 

o 

(AA) 2 . To find a characteristic value for AA, we solve for HWHM (A) 
above: 


(AX) 2 
o 

HWHM (A) = AX 


T 

3.16 x 10 6 
T 

3.16 x 10 6 


For T = 1000 K, HWHM (A) = 0.01776 


x 10 
from 
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